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It is proposed that the upper pseudogap phase (UPP) observed in the high-T c cuprates correspond 
to the formation of spin singlet pairing under the bosonic resonating-valence-bond (RVB) descrip- 
tion. We present a series of evidence in support of such a scenario based on the calculated magnetic 
properties including uniform spin susceptibility, spin-lattice relaxation and spin-echo decay rates, 
which consistently show that strong spin correlations start to develop upon entering the UPP, being 
enhanced around the momentum (it, it) while suppressed around (0, 0). The phase diagram in 
the parameter space of doping concentration, temperature, and external magnetic field, is obtained 
based on the the bosonic RVB theory. In particular, the competition between the Zeeman splitting 
and singlet pairing determines a simple relation between the "critical" magnetic field, -ffpc and 
characteristic temperature scale, To, of the UPP. We also discuss the magnetic behavior in the lower 
pseudogap phase at a temperature T v lower than To, which is characterized by the formation of 
Cooper pair amplitude where the low-lying spin fluctuations get suppressed at both (0, 0) and (ir, 
7r). Properties of the UPP involving charge channels will be also briefly discussed. 

I. INTRODUCTION 

The pseudogap phase has been widely regarded as an essential integral part of the cuprate superconductors with 
extensive experimental support i The underlying physics of such a phase has been a central focus in the study of high- 
Tc problem for many years, and yet no consensus has been reached concerning its nature due to the very complexity 
of pseudogap phenomena observed in magnetic, transport, single-particle, and optical channels. More and more 
experimental evidence in recent years further indicates the existence of two kinds of pseudogap regimes at different 
temperatures 

Among various theoretical proposals for the pseudogap phase, the RVB idea£ is uniquely interesting, which actually 
"predicted " 78 9 10 the existence of a pseudogap state in doped Mott insulators before experiment. In the RVB picture, 
neutral spins form the singlet pairs that are condensed as a spin liquid. The density of states of spin excitations for 
such a system generally gets suppressed at low energy, exhibiting a pseudogap feature as it costs energy to break up 
the RVB pairs to create spin excitations. In such a scenario, there is usually no gap in the pure charge degrees of 
freedom, and the pseudogap phenomena observed in the charge transport, angle resolved photoemission spectroscopy 
(ARPES), and tunnelling experiments, are all indirectly attributed to the appearance of the spin gap in the spin 
degrees of freedom^ For instance, in the charge transport the strong scattering between the charge carriers and 
low-lying spin fluctuations becomes weakened because of the reduction of the latter in the pseudogap regime; The 
pseudogap feature exhibited in the ARPES and tunnelling measurements may be also interpreted as due to the opening 
of a pseudogap associated with the spin degrees of freedom. 

However, the original RVB description, known as the fermionic RVB (f-RVB) since it involves the pairing of fermionic 
"spinons" (neutral 5 = 1/2 object) j£ also suffers some notable inconsistency with the experiment. Note that the 
pseudogap phenomenon has been found in the underdoped regime of the cupratesi where the antiferromagnetic (AF) 
correlations are usually quite strong. But in an f-RVB description, the AF correlations remain intrinsically weak, 
even at half-filling, where the AF long range order (AFLRO) develops in the cuprates at low temperature. Here 
the key issue is not about whether one can construct an AFLRO in the RVB background, which may be easy to 
incorporate by a simple mean-field order parameter. 11-12 But the crucial and highly nontrivial issue is whether the 
• *h , whole low-lying AF spin correlations are intrinsically and sufficiently strong in an RVB stated, and whether they are 
capable of continuously growing with reducing temperature or doping as has been clearly manifested experimentally 
in, say, the NMR spin-lattice relaxation rates (see the analysis in Refsi^i£*i£iiL 1 £). 

Since the pseudogap phase, which involves high-energy/temperature and short distance physics, may be properly 
considered as an unstable fixed point stated with intrinsic instabilities towards AFLRO or d-wave superconductivity 
at low temperatures, the importance of its correct description, with regard to the latter, is like that of a Fermi 
liquid with regard to the BCS superconducting state. In other words, finding an accurate and correct description 
of the pseudogap phase will be rather important for constructing a sensible low-energy theory for describing the 
low-temperature AF and superconducting phases in the cuprates. 

In this paper, we show that there does exist a desirable candidate for characterizing the pseudogap phase based 
on the RVB picture, in which strong AF correlations are present as an intrinsic and predominant feature. Such an 
RVB state, known as the bosonic RVB (b-RVB) state, 19 differs from the usual f-RVB states by that it works very 
well at half-filling in describing the AF correlations over a wide range of temperature, including zero temperature 
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FIG. 1: The global phase diagram in the b-RVB description (Refill). The upper pseudogap phase (UPP) is characterized by the 
bosonic RVB pairing order parameter A s at T < To whose properties are the main focus of this paper. The antiferromagnetic 
ordered phase (AF), the lower pseudogap phase at T < T v , and the superconducting phase (SC) at T < T c all happen on top 
of this UPP at low doping. 

where an AFLRO naturally emerges. Fig. schematically illustrates the global phase diagram for such a b-RVB 
theory, 20 where an upper pseudogap phase (UPP) below the characteristic temperature Tq is characterized by the 
formation of singlet pairing of neutral spins as denoted by the b-RVB order parameter A s . As illustrated by Fig. 
1, the low-temperature phases, including a lower pseudogap phase below T v (also known as the spontaneous vortex 
phase 21 ), a d-wave superconducting phase at T c < T v beyond a critical doping concentration x c , and an AFLRO 
phase near half-filling, can all be regarded as the results of the low-temperature instabilities from such an UPPiS . 

The main focus of the present work will be the nature of the UPP itself. We will examine the detailed behavior of 
uniform spin susceptibility, spin-lattice relaxation rates, and spin-echo decay rate based on the b-RVB theory, which 
will reveal that, as one enters the UPP from above To, the spin correlations change qualitatively. At the mean- field 
level, localized spins are essentially uncorrelated at T > To, resulting in a Curie- Weiss-like behavior in the uniform 
spin susceptibility. Below To, however, finite-range spin correlations start to develop, predominantly around the AF 
momentum Qaf = (""ji"); with the weight being transferred from the momentum Qo = (0,0). The latter effect leads 
to the reduction of uniform susceptibility at T < Tq. We show that such magnetic behavior in the UPP is quite 
consistent with the experimental measurements in the cuprates. 

The above results clearly indicate that the UPP corresponds to a crossover from a weakly correlated localized spin 
system at higher temperature into a strongly AF correlated spin liquid at lower temperature. This picture is thus in 
sharp contrast to the f-RVB picture of the pseudogap, where the low-energy spin excitations, either around Qaf or 
Qo, all get suppressed with the opening of the pseudogap. This latter behavior rather resembles the lower pseudogap 
phase of the b-RVB theory at T < T v (in Refi2£ it is simply called the pseudogap phase) than the UPP. But even the 
distinction between the pseudogap state of the f-RVB theory and such a lower pseudogap state of the b-RVB theory 
is very significant as has been discussed in Ref^: the former is exchange energy driven while the latter is kinetic 
energy driven. Here the lower pseudogap phase corresponds to the formation of Cooper pair amplitude but is short 
of superconducting phase coherence, which exists between T v and T c and can be regarded as a vortex liquid stated 
It was previously pointed out in Ref^i that such a spontaneous vortex phase should coincide with the experimentally 
discovered Nernst region 2 ^ in the high-T c cuprates. 

The quantitative phase diagram of the UPP is determined in the three-dimensional parameter space of temperature, 
doping concentration, and external magnetic field. The latter introduces the competition between the Zeeman spin 
splitting and singlet spin pairing. The mean-field theory will predict a simple proportional relation between the 
"critical" magnetic field Hpq, at which the UPP is destroyed, and Tq in zero field. A comparison with experiment 
will be made. 

The remainder of the paper is arranged as follows. In Sec. II A, the bosonic RVB description is briefly reviewed. 
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In Sec. II B, the definition of the upper pseudogap phase is given and its phase diagram is determined. In Sec. II C 
and D, its magnetic properties are systematically investigated in the framework of the b-RVB theory. In Sec. II E, 
we further briefly discuss the lower pseudogap phase and related magnetic behavior. Finally, Sec. Ill is devoted to 
conclusion and discussion, where we also briefly discuss the qualitative behavior of charge channels in the pseudogap 
phase within the bosonic RVB description. 



II. UPPER PSEUDOGAP PHASE IN THE BOSONIC RVB THEORY 



A. Bosonic RVB description 

The so-called b-RVB stated is underpinned by a bosonic RVB order parameter A| ■ over a wide range of temperature 
and doping as schematically shown in Fig. 1, which defines the UPP. This UPP (before the emergence of low- 
temperature AF and superconducting instabilities) will be the main subject to be examined in the present work. In 
the following we shall first discuss its "mean-field" description based on the t — J model. 

In the phase string representation 2 ^ of the t — J model (see Appendix A) , a natural "mean-field" decoupling of the 
superexchange term Hj is given as follows*^ 

hj^h s = ~j2 ^ iaA HA-° + H - c - + const - (!) 

where the b-RVB order parameter Af^ is defined in terms of the bosonic spinon annihilation operator bi a by 

A ^=E( e "^^-)- (2) 

cr 

At half-filling, Af ■ is equivalent to the usual Schwinger-boson mean-field order parameter as Aj-j = 0, with Eq.QJ 
reducing to the Schwinger-boson mean-field Hamiltonian 25 which describes the AF correlations fairly well in the regime 
of Afj ^ at T < To = 0.91J/k B . Away from half-filling, a topological link field emerges in the Hamiltonian (JTJ 
as well as in Eq. , which represents the influence of the nonlocal phase string effect induced by the hole hopping 24 
on the spin degrees of freedom. It is related to the hole density by the following gauge invariant relation^ 

E 4 = -E< (3) 

<ij>ec ies c 

where C denotes an arbitrary loop on the square lattice that encloses a region £c, and is the holon number 
operator at site 

The "mean- field" Hamiltonian Q is by nature a gauge model, which is invariant under the U{1) transformation: 
bier —> bi a e la ' e * and A^ — > A^ + (9i — 6j). It can be easily shown that the spin rotational symmetry is respected 
by Eq. by verifying [H s , S] = where the definition of the spin operator S in the phase string representation is 
given in Appendix A. In the low-temperature regime, where the bosonic holons in the b-RVB theory will experience 
the Bose condensation such that one may approximately treat A^ as describing a uniform flux of the strength Sir 
per plaquette (S denotes the doping concentration of holes) In the high-temperature regime, A^ may be treated 
as describing randomly distributed static ir flux tubes of concentration S since the holons will behave like incoherent 
objects there^ Different from a usual Jordan- Wigner phase, the gauge field A^., which is seen by spinons, is attached 
to an independent degree of freedom, holons, and therefore the above approximations are reasonable. 

In both limits, the bilinear form of Eq.Q in terms of the bosonic spinon operator bier can be diagonalized by the 
following Bogoliubov transformatiomiS 

h 'io = ^["mj(')7m(i ~ iw(i)7m-<r] ( 4 ) 



with 



v mc r(i) = v m w m<y (i) (5) 



where w ma (i) satisfies the following eigen equation 

" 2 



£ ro uw(i) = -~ E A^e-^uwtf) (6) 

j=nn(i) 
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in which j = nnli) denotes the four nearest neighbors (nn) of site i. One has u m — ^Tfyir — l~ 1 an d v m — 



sgn(£ m )^y-gr 1, where E m = ^/X 2 — £ m * s the s pi non spectrum. The Lagrangian multiplier A is determined 

by enforcing the average constraint (J2<j "iJ^ic) — 1 — which leads to 



m 

in) 

where j3 = l/ksT and the last equation is obtained by the self-consistent condition (J2J) for the RVB order parameter. 
Note that a Bose condensed term n b BC related to the AFLROiS at T — and half- filling has been dropped in Eq. 1131) 
since we shall mainly be interested in the high-temperature behavior. 

All the nontrivial effect of doping is reflected in the eigen equation where the phase string effect induced by 
hopping enters via the topological gauge field A^.. Note that the spinon wave function w ma (i) should vanish at the 
hole sites (where A|^ = 0) due to the no double occupancy condition. Previously^ such an equation has been solved 
by a simple mean-field choice A*- [j — nn(i)]= A s with the relaxed constraint condition such that the spinons can 
go any sites. In the following we still relax the no double occupancy constraint on average in Eq.©, but with a 
replacement of 



J -» Je ff = (1 - 2gS)J (9) 

to represent the average effect of the reduction of the superexchange coupling around holes: if holes are static, each 
of which will simply break two nn links in each direction such that g = 1 in the dilute hole limit (when the average 
hole-hole distance is much larger than the nn links). Generally g > 1 for a moving hole since the suppression of Af^ 
around each hole extends more than the four broken nn bonds (e.g., from the singular twist by around each hole). 
One may thus approximately rewrite Eq.© as 

£mWm*(i) = -J s e ~ iCrA ^ w ma{j) (10) 

j=nn(i) 



in which 



and Eq.JHJ) can be consistently rewritten as 



JcffA s 



A s = — Y Jl^coth^. (11) 



Note that the same mean-field equations have been obtained in Ref^ with a slightly different definition, i.e., with 
A s replaced by A s /(1 - 2gS) = A{ (with g = 1) in Ref-iS. 



B. Upper pseudogap phase 

The UPP is defined by the formation of the b-RVB pairing with A s ^ 0. Its high-temperature boundary at A s = 
is depicted by a characteristic temperature To as illustrated in Fig. ^ In the following, we determine it based on the 
mean-field theory outlined above. 

According to Eq. Q , one finds 

2-* = coth-A_ (12) 
2fc B T 

by noting that £ m — > at A s — > 0. Consequently, Eq.JTT} reduces to 1 = (2 — S)Jjj ^ m (£ m /A s ) 2 / J e ffA, which gives 
rise to 

2 5 

A = Jeff (13) 
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by further identifying 

l]£(£ m /A s ) 2 = J 2 ff (14) 



N 

in terms of Eq. I|1(J|I . Finally one obtains 




k BT = Jeff ■ (15) 



It is interesting to note that the gauge field does not explicitly appear in T because of the general relation 
(fH|) . indicating that the phase string effect, which is crucial to the low temperature (low energy) physics, actually 
plays no role in determining the phase boundary of the UPP. Of course, the hopping effect still enters in To via the 
renormalized factor g in J c fr which shall be the only adjustable parameter depending on the detailed physics of local 
hopping. We find that the phase diagram and magnetic properties to be studied below are actually not sensitive to 
g except for the characteristic concentration xrvb at which J e g is extrapolated to zero at T = 0. 

Fig. |21 shows To (solid curve) as a function of 8/xbyb with J = 1350 K. The experimental data determined by the 
uniform spin susceptibility measurement in LSC0 26,2 £ (see the discussion in the next section) are shown by the full 
squares, which are in good agreement with the theoretical curve. Furthermore, the open circles are independently 
determined from the c-axis transport in the overdoped regime (see discussion below). Note that the theoretical 
curve T) versus 5/xrvb in Fig. 2 is not sensitive to the choice of g. Furthermore, one can use the above experimental 
dat a 26 i 27 i 28 to fix xbvb at £rvb = l/2g = 0.25 (g = 2). We shall then choose itrvb = 0.25 throughout the rest of 
paper without any more adjustable parameter. 

For the b-RVB origin of the UPP, the Zeeman splitting due to the external magnetic field can effectively destroy 
the singlet pairing of spins in the strong field limit. Since the orbit part of the neutral spins does not couple to 
the external field directly, the Zeeman splitting will be the only direct field effect on the RVB background. It thus 
provides a direct probe of the RVB nature of the UPP in, say, the overdoped regime, where the critical field strength 
may be within the experimental accessible range. In the following we consider the Zeeman effect in the UPP. 

Apply an external magnetic field H along a spin z-axis (which is not necessarily perpendicular to the lattice plane). 
A spin Zeeman energy term should be then added to H s in Eq.Q: 

- 2/i B J2 S * H = -VBHj^^rnalma (16) 

i a 

Consequently the spinon excitation spectrum is modified by 

E° m = E m - <jii b H (17) 
which now explicitly depends on the spin index a. Then the mean-field equation i|12f> at A s — > is modified to: 

2 - S = - Vcoth-^— (18) 
2 ^ 2Tk B 

with E a = A — (TfiBH, while Ea. ()l^j) remains the same. From these equations, we can easily obtain the following 
relation between To at zero field and the zero-temperature "critical" field H PG at which A s vanishes: 

with the coefficient only weakly dependent on the doping concentration. In Fig. [21 M-B-^pg/^ 5 In 3 is plotted as the 
dashed curve which scales with the zero-field To fairly well, which predicts 

(i B H% G ~l.lk B T . (20) 

In general, the temperature dependence of the "critical" field -ffpc(T) can be obtained based on Eqs. (|18fl and l|13f) . 
In the inset of Fig. [21 -ffpc versus T at S = 0.22 is shown together with the experimental data obtained from the 
c-axis magneto transport measurements f 2 ^ We see that the high-temperature experimental data (open circles) fit the 
theoretical curve very well without any additional adjustable parameter. Furthermore the zero-field Tq determined by 
the same experiments is also in good agreement with the theory as shown (open circles) in the main panel of Fig. [21 
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FIG. 2: The characteristic temperature To of the UPP versus Solid line: the present theory; Full squares: determined 

from the uniform spin susceptibility \s in LSCO compound;— Open circles: determined from the c-axis magneto-resistivity 
(p c ) measurement in Bi-2212 compound?*- The dashed line shows the scaling relation of the zero-temperature critical field Hp G 
with To as predicted by the theory. Inset: the critical field -Hpg as a function of temperature at S = 0.22. The experiment data 
from the c-axis transport in Bi-2212 (Refi»&) are also shown by the open and full squares. 

as mentioned above. One may also note that the experimental Hpq(T) starts to deviate from the theoretical curve in 
the inset (full squares) as the temperature is further lowered and saturated to approximately the half of the predicted 
number (which implies [IbH^q ~ fcsTo/2). However, we would like to point out that such a deviation occurs only for 
those data (full squares) which have been obtained by extrapolation in the experimental measurements^ and therefore 
may not be as reliable as the higher temperature ones (open squares) in the inset of Fig. We believe that further 
experiments be needed in order to convincingly verify (falsify) the present prediction l|2(J|) . 

Finally the boundaries of the UPP in the three dimensional parameter space of doping concentration, temperature, 
and magnetic field determined based on the b-RVB mean-field theory are shown in Fig. [3J 



The RVB nature of the UPP is clearly manifested in the magnetic properties. We first consider uniform spin 
susceptibility \s in the following, which can be easily derived based on the above mean-field description in the 
presence of magnetic field. In terms of the total spin magnetic moment 



[where n(u>) — l/(e^ u — 1) is the Bose distribution function], the uniform spin susceptibility \s per cite as defined by 
^ s = Wh l#^o is found by 



C. Uniform spin susceptibility 




(21) 



Til 




(22) 



n i 



Since we are mainly interested in the high-temperature regime well above the superconducting phase, the holes are 
incoherent objects such that can be approximately treated as describing randomly distributed tt flux tubes of 
concentration S as discussed before. Then we can numerically calculate Xs based on the mean-field equations given 
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FIG. 4: (a) The calculated uniform spin susceptibility Xs scaled with the maximum x™ ax at To versus T/To, which follows an 
approximately doping-dependent curve. Inset: The experimental data in Refi— *— which collapse into a universal scaling curve 
plotted in the same fashion as in the main panel, (b) The theoretical Xs at half-filling (solid) and the one obtained by the high 
temperature series expansion (HTSE). The latter fits the experimental scaling curve in the inset of (a) very welli— ■ *— 



in Sec. II A, which is averaged under different random configurations of A^. A similar computation has been done 
before to explore the crossover from the lower to upper pseudogap phases^ but not in the region up close to To. 
(In the following calculations, the largest sample size is 32 x 32 lattice and the sample size is not very important as 
mainly the high temperature properties are concerned.) 

The calculated Xs is presented in the main panel of Fig. 0] (a) at different doping concentrations. Note that Xs 
reaches a maximum value xT ax a ^ temperature Tq where the RVB order parameter A s vanishes. At T > Tq, Xs follows 
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a Curie-l/T behavior as spins become free moments at the mean-field level. The curves in Fig. 01(a) are presented as 
Xs/x™ ax versus T/T Q , which approximately collapse onto a single curve independent of doping. For comparison, the 
inset shows the experimental data obtained in LSCO compounds which are plotted in the same way as in the main 
panel with a very good collapsing onto a universal scaling curve; 26,27 And the peak positions decides the experimental 
pseudogap temperature To's at different dopings, which are presented in Fig. [2] 

In Fig. 0](b), the calculated %s versus T at 6 = is shown together with the high temperature series expansion 
(HTSE) result (note that here the calculated \s is rescaled by a 2/3 numerical factor as used in the Schwinger-boson 
approach to restore the sum rule2£). It is noted that the experimental scaling curve actually coincides with the half- 
filling HTSE very welli^ 6 ^ 7 . Thus one can clearly see the overall qualitative agreement between the bosonic RVB theory 
and the experiment from Figs. |H(a) and (b). Note that the mean-field Xs deviates from the HTSE result prominently 
around T where the latter is a much smoother function of T. It reflects the fact that T is only a crossover temperature 
and the vanishing A s does not represent a true phase transition. Obviously, the amplitude fluctuations beyond the 
mean-field A s have to be considered in order to better describe Xs in this regime. To determined in the mean-field 
theory is quite close to the HTSE result, indicating the crossover temperature itself can still be reasonably decided 
by the mean-field bosonic RVB description given above. The comparison of To between the theory and experiment 
has been already presented in Fig. [21 and discussed in the previous section. 



D. Spin-lattice relaxation rate and spin-echo decay rate 

The NMR spin-lattice relaxation rate of nuclear spins is determined by 3Q 

1 2k B T \ - 2 

where the form factor T(q) 2 comes from the hyperfine coupling between nuclear spin and spin fluctuations. For 
example, for planar 63 Cu(2) nuclear spins in the cuprates, with the applied field perpendicular to the Cu02 plane, 
the form factor T(q) 2 is found to be*i£*i& 

63 F(q) 2 = [A ± + 2B(cos q x a + cos q y a)} 2 (24) 

where the hyperfine couplings A± and B are estimated as A±/B ~ 0.84, B ~ 3.8 x 10 _4 meV— [These coefficients 
may slightly vary among YBCO and LSCO compounds]. For planar 17 0(2) nuclear spins, one has 

17 F(q) 2 = 2C 2 [1 + -(cosg x a + cosmos)] (25) 

with C ~ 0.87-B. Due to the fact that 17 T(q) 2 vanishes at the AF wave vector Qaf = (^,^), while 63 F(q) 2 is 
peaked at Qaf, a combined measurement of l/ 63 Ti and l/ 17 Ti will thus provide important information about the 
AF correlations at low frequency lu — > 0. 

Based on the bosonic RVB mean-field equations outlined in Sec. II A, the spin-lattice relaxation rates, l/ 63 Ti and 
l/ 17 Ti, for the planar copper and oxygen nuclear spins, can be straightforwardly computed as shown in Appendix 
B. The results are presented in Fig. It shows that the ratio 17 Ti/ 63 Ti, which is a constant above To, starts to 
increase with reducing temperature below To. At lower temperature, T/Tq < 0.5, such a ratio arises sharply. For 
example, 17 Ti/ 63 Ti diverges at S = as a true AFLRO exists at T = 0; And it can still reaches about 100 in the 
low temperature limit at 5 — 0.125, all qualitatively consistent with the experimental observation in the cupratesiii 
As pointed out above, such behavior clearly demonstrates that strong low-lying AF correlations around Qaf develop 
in the UPP, leading to the simultaneous enhancement of l/ 63 Ti and the cancellation in l/ 17 Ti. In the inset of Fig. 

i/w^j* j-^g k een plotted, which is also qualitatively consistent with the experiment jiLiLiS but deviates from the 
conventional Korringa behavior l/ 63 TiT ~ const for a Fermi liquid system, thanks to the strong AF fluctuations 
of spins in the UPP below To^i By contrast, the ratio 17 Ti/ 63 Ti in an f-RVB mean-field state (the it flux phase) 
at half-filling remains flat over the whole temperature region as shown by the dashed line in Fig. indicating the 
absence of any significant AF correlations around Qaf in its pseudogap regime. In combination with the reduced 
uniform spin susceptibility— one sees that in the pseudogap of the f-RVB, the low-energy spin excitations, cither 
around Qaf or Qo = (0, 0), all get suppressed with the opening of the pseudogap. 

The above results clearly show that the UPP in the bosonic RVB state corresponds to the crossover from a weakly 
correlated localized spin assembly at higher temperature into a strongly AF correlated spin liquid at lower temperature. 
The peculiar feature of the bosonic RVB description is that although the formation of bosonic RVB singlet pairing 
suppresses the spin correlations at Qo below To, it also leads to the enhancement of the low-energy spin correlations 
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FIG. 5: 17 Ti/ 63 Ti vs. temperature at different doping concentrations in the upper pseudogap phase of the b-RVB state. The 
dashed fine shows the resuit of an f-RVB state (ir flux phase) at half-filling. The inset shows the non-Korringa behavior of 
l/ 63 TiT in the b-RVB state at various dopings (the symbols are the same as in the main panel). 



near AF momentum Qaf- Such a feature of the UPP is significantly different from the pseudogap phase in the f-RVB 
mean-field description, but is strongly supported by the NMR measurements ^iiSiiSiiLiSi 

Finally, we further examine the spin-echo decay rate 1/T 2 g, which is related to the static AF correlations via the 
real part of spin susceptibility function by>22 



1 

Tig 



0.69 l_ 

8H 2 {fllef 



l^F e// (q) 4 X L(q) 2 



^E^//(q) 2 xL(q) 



(26) 



where the factor F e g(q) is 



^eff(q) = A\\ + 2£?(cos q x a + cosq y a) 



(27) 



with Ay ~ — AB, such that F e gp(q) is peaked at Qaf and vanish at Qo- 

Similar to 1/Ti, the detailed expression of I/X2G m the bosonic RVB theory is given in Appendix C. In FigEl the 
calculated I/T2G a t different doping concentrations show that I/T^g begins to increase with reducing temperature 
below To. Such behavior has been also observed in the experiment j&S&iSk which once again clearly illustrates the 
picture that the strong AF correlations start to develop in the UPP. 



E. Lower Pseudogap Phase 



So far we have been focused on the UPP, which is the high-temperature phase in the bosonic RVB description. As 
stressed in the Introduction, there can be several different low-temperature phases growing out of this background (Fig. 
\)M> One particular phase we wish to discuss below is the so-called spontaneous vortex phased which can be properly 
classified as the lower pseudogap phase in the present approach. In this phase, the holon condensation occurs and 
Cooper pair amplitude forms, but the system is still short of superconducting phase coherence which can be regarded 
as a vortex liquid state due to the presence of unpaired spinon-vortex composites. 21 These spinon-vortices contribute 
to the Nernst effect and therefore the lower pseudogap phase should coincide with the Nernst region discovered 2 ^ 
experimentally in the cuprates. Recently the electromagnetic response of such a vortex liquid phase has been also 
discussed 36 based on a different RVB approach. 

Previously the magnetic properties in this phase has been discussed 2 ^ in the context of exploring the driving 
mechanism in comparison with the pseudogap phase in the f-RVB state. In the following, we focus on the magnetic 
behavior of this lower pseudogap phase and make contrast with that of the UPP discussed above. 
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FIG. 6: I/T2G vs. temperature in the upper pseudogap phase below To at different doping concentrations. 




FIG. 7: Uniform spin susceptibility in the lower pseudogap phase at different dopings including half-filling. The left insert 
shows l/ B3 Ti and the right insert I/T2G with the same symbols as in the main panel. 



The mean-field equations are the same as in the UPP, but due to the holon condensation the gauge field in 
Eq. H10JI can be treated as a uniform flux of strength 7r<5 per plaquette 2 ^ as discussed in Sec. II A. In the main panel 
of Fig. [7] the uniform spin susceptibility shows a true "spin gap" behavior, in contrast to the "scaling" curve shown 
in the UPP in Fig. 4 where %s in the doped regime roughly behaves like that at half-filling — in the latter case Xs 
saturates to a constant at T — 0. In the lower pseudogap phase, these Xs's can drop below that at S = and vanish 
at T = as shown in Fig. Such a lower pseudogap behavior has been indeed observed experimentally. 2,9 

Furthermore, in this lower pseudogap phase, l/ 63 Ti also decreases with temperature (see the left inset of Fig. 7), 
as opposed to the behavior in the UPP, indicating the appearance of the spin gap over whole momenta. On the other 
hand, although the low-energy spin fluctuations are gapped, the static AF spin-spin correlations as described by the 
real part of spin susceptibility function still remain, as reflected by I/T2G shown in the right inset of Fig. 7, where 
the monotonic increase of I/T2G in the UPP (Fig. 6) is replaced by the saturation at lower pseudogap phase. This 
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feature has also been observed experimentallyiSiMi^ 

III. CONCLUSION AND DISCUSSION 

In this paper, we have systematically analyzed the magnetic characterizations of a high-temperature intrinsic 
phase of the bosonic RVB state, which is described by the formation of the bosonic RVB order parameter at a 
temperature below the characteristic To, but still higher than those for low-temperature orders, including AFLRO 
and superconductivity, to emerge. 

Such a phase exhibits the pseudogap features that match those of the upper pseudogap phase in the high-T c 
cuprate superconductors very well. The key feature in the crossover to the UPP from above To is the onset of the 
development of strong AF spin-spin correlations, which remain rather weak at T > To where the system resembles 
more an ensemble of uncorrelated localized spins. This explains why experimentally the uniform spin susceptibility 
shows an approximately Curie- Weiss behavior at T > To, reaches a peak at To, and then gets reduces below To as 
the weight of the spin-spin correlations at the momentum (0,0) being transferred to (tt, it), in contrast to an equal 
weight distribution above Tq. It further explains why the spin-lattice relaxation rate gets enhanced for the planar 
copper nuclear spins whereas reduced for the planar oxygen nuclear spins below To, and why the spin-echo decay 
rate increases with the decreasing temperature; Clearly the development of the AF correlations is the underlying 
mechanism here. 

We emphasize that the formation of spin singlet pairing and the onset of AF correlations at the same time are 
not always true. In the f-RVB description of the slave-boson approach, the formation of the f-RVB order parameter 
actually leads to the reduction of the spin-lattice relaxation rates for both the copper and oxygen nuclear spins at low 
temperatures. This is because the spin pseudogap opens for both the ferromagnetic and AF correlations. By contrast, 
this case occurs in the b-RVB description only at a lower temperature when the system enters the lower pseudogap 
regime, characterized by the formation of Cooper pair amplitude in the so-called spontaneous vortex phaseSi which is 
a vortex liquid state, short of superconducting phase coherence. A comparative study of this lower pseudogap phase in 
the b-RVB theory and the pseudogap phase in the f-RVB theory has been given in Ref^, where two opposite driving 
mechanisms: kinetic vs. superexchange energy driven, have been identified. 

The phase diagram of the UPP has been determined by a generalized "mean-field" description in the b-RVB theory 
in the parameter space of temperature, doping, and magnetic field. Since the UPP essentially reflects the spin singlet 
pairing, the Zeeman splitting competes directly, which results in a quantitative prediction for experiment as discussed 
in Sec. II B. 

A central consequence of an RVB (spin singlet) description of the pseudogap phase is that the pseudogap represents 
the suppression of the low- lying spectral weight in spin excitations, but not in charge excitations, as pointed out in 
the Introduction. In the following, we briefly outline the scenario about what happens to the charge channel when 
one enters the UPP in the bosonic RVB theory. We shall leave the more quantitative investigation in future work. 

In order to see how the charge transport is affected by the spin fluctuations, we first note that in the b-RVB theory, 
the charge degrees of freedom (holons) are described by the following effective HamiltoniantiS 

H h = -t h e iA vh\hj + h.c. (28) 

<ij>cr 

where the bosonic holons, created by h\, interact with the gauge field Afj associated with the spin degrees of freedom. 
Similar to the definition of the gauge field A 1 }- in Eq. J3J) , A\^ , which is introduced in the phase string representation 
(see Appendix A), satisfies 

E ^ = 7r E {<-<) ( 29 ) 

<ij>ec zgSc 

where C is an arbitrary path. Physically Afj describes ±7r flux tubes bound to f(j.) spinons as seen by holons. So at 
T > T), uncorrelated localized spins imply a maximum scattering to the holons according to Eas. (|28|) and (|29[1 . By 
forming the RVB pairing below To, one can easily understand that the fluctuations in Afj will be effectively reduced, 
and so does the scattering to the holons according to Eqs.J2HJ an d (|29|l . leading to a pseudogap feature in the charge 
transport without involving a charge gap. Note that when the temperature is further reduced to T v , where the holons 
gain the phase coherence and become Bose condensed, the system enters the lower pseudogap phase (spontaneous 
vortex phase) in which the effect of the holon condensation will feedback to the spinon part via A^ in Eq. and cause 
the lower pseudogap phenomenon in the spin part as discussed in Sec. II E. Finally, the quasiparticle excitation can 
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be regarded as a recombination of holon, spinon, and phase string, in the superconducting phase. 3 It has been argued 
that the deconfincmcnt occurs above T c , and the composite structure is expected to shows up in both the lower and 
upper pseudogap phases, and the pseudogap feature is thus believed to be associated with that in the spinon degrees 
of freedom. 
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APPENDIX A: PHASE STRING FORMULATION 



The t - J model 



t (4^ + h.c.) + jJ2(Si- Sj - -mnj) (Al) 



H t -j 

<ij><J <ij> 



may be reformulated by using the phase string decompositio: 



.24 



c ia = {~u) l h\b ia e<" na (A2) 

where hi and bi a are all bosonic fields. Here Q^™ 9 ig a non-local phase factor to restore the fermionic statistics 
of the electron operator, and can be expressed as Q^*"™ 9 = — a^] with $| = X^i^MO (Ea ra L — l)> 

— J2i^i #i(O n f - Here 9i(l) is defined as an angle 9i(l) — Im ln(zi — z\) with = X{ + iyi representing the complex 
coordinate of a lattice site i. The resulting Hamiltonian H t -,j — H t + Hj reads 

<ij>a 

Hj = J2 (e^blbt^'^b^b^ (A3) 

under the no-doublc-occupancy constraint 

h\hi + Y,blb ir , = l . (A4) 



In the new Hamiltonians, <^ is a 7r flux link variable, while Afj and are constrained by the following conditions: 

c lec 

£4 = nj^n? (A6) 



where C is an arbitrary counterclockwise closed path. Finally, in the phase string representation, the spin operators 
can be easily reexpressed according to the decomposition (|A2(1 as 



S- = (-mlk-^ . (A7) 
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APPENDIX B: 1/Ti FORMULATION IN THE BOSONIC RVB REPRESENTATION 

1/Ti defined in (|23|) can be reexpressed in terms of the real-space spin correlation function as follows 

1 _ 2k B T ^ x'Lihj^) 



uJ 



(Bl) 



where My is the Fourier transformation of F(q) in real space: 

63 



Mij = (Ai+4B 2 )S lJ +2A ± BY,^+f l +B 2 ]T <W )+ ,y (B2) 
where f] = ±x, ±y and 

17 M^ ee 2C 2 + ~ ]T tfij+fl) . (B3) 

From Mij we see that only up to the next-nearest-neighbor spin correlations are involved in the spin-lattice relaxation 
rates of 63 Cu and 17 O nuclear spins. 

In the Bosonic RVB mean- filed theory, by using Eq.(0} the dynamic spin susceptibility can be expressed as 



UJ 

where 



= Gj. + (-iy-iG+ (B4) 



G% = ^E' K ^(iJ) f-^ 1 ) (pL,) 2 S(E m - E m ,) (B5) 



2^ mm'V^M Q En 
mm' 

in which denote the summation of m with £ m > and 

K mmf(hj) = ^w*^(i)w mrT {j)w* m , a (j)w m ^{i) (B6) 



By noting that p mm , = l,p^ m / = X/E m at B TO = E m i, we can further reexpress 1/Ti in the following form 



D m + §T D r. 



T x 3g^ B N< 

where the density of states p(E m ) — (2/N)Y^ m S(E m — -E m <) and the coefficient, D^, is defined by 



(B7) 



n ± Em d mm ,6(E m E m i) 
Z^m d <A» _ E m >) 

with 

^ = f N E (i,i)(Tr^M y (B9) 

In Ea. llB7() a numerical factor 2/3 is also added just like the uniform spin susceptibility as noted in the main text as 
at half-filling. 25 The final result will be an average over different random configurations of A^- due to the incoherent 
distribution of holes in the UPP (In the lower pseudogap phase, by contrast, the holon condensation leads to a uniform 
flux distribution of A^ and no such average is needed). The calculation is done on a 32 x 32 lattice and results are 
presented in Fig. [5J 31 
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APPENDIX C: SPIN-ECHO RELAXATION RATE 

The spin-echo relaxation rate I/T2G is defined in Ea. (|26() . In the b-RVB theory, the real part of the static suscep- 
tibility, x'zz (q) i can be expressed as 



x' z M = x' z + M + x' z -A<D (ci) 



with 



x2(q) = |x^V- m ,(q) 



± , 2 (n(E m >) - n{E m )) ± 2 (1 + n(E m ) + n{E m ,)) 



\"mm / tp tp 1 V nun J 



Tp fp ' \ mm / Tp I tp 



(C2) 



where 



Pmm' — ^m^m' i ^m^m' 



mm 

^mm 1 ~ u m v m' i V m U m i (C3) 

and 



N 



The numerical calculation is similar to that for 1 /Ti . 
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